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ABSTRACT 

The  motion  of  a  two-dimensional  bubble  rising  at  a  constant  velocity  In 
an  unbounded  fluid  Is  solved  by  series  truncation*  It  Is  assumed  there  Is  a 
wake  of  stagnant  liquid  extending  to  Infinity  below  the  bubble.  Both  the 
effects  of  gravity  g  and  surface  tension  T  are  taken  Into  account.  It  Is 
shown  that  the  problem  Is  characterized  by  a  continuum  of  solutions  for  T  =  0 
and  by  a  discrete  set  of  solutions  when  T  >  0.  In  addition  a  unique  solution 
Is  obtained  In  the  limit  as  the  surface  tension  approaches  zero.  The 
corresponding  profile  of  the  bubble  Is  found  to  be  In  good  agreement  with 
experimental  data. 
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SIGNIFICANCE  AND  EXPLANATION 


Recent  numerical  computations  on  nonlinear  free  surface  flow  problems 
(Vanden-Broeck^'^'^)  have  uncovered  an  unexpected  effect  of  surface  tension. 

It  has  been  found  that  some  problems  are  characterized  by  a  continuum  of 
solutions  when  surface  tension  Is  neglected  and  by  a  discrete  set  of  solutions 
when  surface  tension  Is  taken  Into  account. 

In  the  present  paper  tra  describe  another^  example  of  such  flows.  We 
consider  the  motion  of  a  two-dimensional  bubble  rising  at  a  constant  velocity 
In  em  unbounded  fluid.  It  Is  assumed  there  Is  a  wake  of  stagnant  liquid 
extending  to  Infinity  below  the  bubble.  Both  the  effects  of  gravity  and 
surface  tension  are  taken  Into  account.  It  Is  found  that  this  problem  Is  also 
characterized  by  a  continuum  of  solutions  when  surface  tension  Is  neglected 
and  by  a  discrete  set  of  solutions  \Aien  surface  tension  Is  taken  Into 
account.  Moreover  It  Is  shown  that  a  unique  solution  Is  obtained  In  the  limit 
as  the  surface  tension  tends  to  zero.  The  corresponding  profile  of  the  bubble 
Is  found  to  be  In  good  agreement  with  experimental  data. 


A  FREE  STREAMLINE  MCX>EL  FOR  A  RISING  BUBBLE 


Jean-Marc  Vanden-Broeck 


1.  Introduction 

We  consider  a  bubble  rising  at  a  constant  velocity  U  In  an  unbounded 
inviscld  fluid.  We  take  a  frame  of  reference  moving  with  the  bubble  and  we 
assume  there  is  a  wake  of  stagnant  liquid  extending  to  infinity  below  the  bubble 
(see  Fig,  1). 

As  we  shall  see  the  problem  is  characterized  by  the  Froude  number 

F  «  0/(gD)''/2  (1) 

and  the  Weber  number 


PO^D 


o  -  .  (2) 

Here  g  is  the  acceleration  of  gravity,  T  the  surface  tension,  p  the 
density  of  the  liquid  and  0  the  width  of  the  bubble  (i.e.  the  distance  between 
the  separation  points  S  and  S'). 

In  this  paper  we  present  numerical  evidence  that  for  (i.e.  T-0) 


there  is  a  solution  for  each  value  of  0  <  F  <  F_  where  F_  ~  0.9.  However, 

c  c 


for  each  value  of  a  ^  there  exists  only  a  countably  infinite  number  of 
solutions.  Each  of  these  solutions  corresponds  to  a  different  value  of  F.  As 
a  tends  to  infinity,  all  these  solutions  approach  a  unique  limiting  solution 
characterized  by  F  =  F*  ~  0.51.  The  corresponding  profile  of  the  bubble  is 
found  to  be  in  good  agreement  with  Collins'^  expejrimental  data. 

.2 


The  present  results  are  similar  to  those  obtained  by  McLean  and  Saffman 
3 


and  Vanden-Broeck for  the  viscous  flow  in  a  Hele-Shaw  cell  and  by  Vanden- 
Broeck^'®'®  for  the  flow  past  a  bubble  in  a  tube.  All  three  problems  are 
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characterized  by  a  continuum  of  solutions  for  T  «  0  and  a  discrete  set  of 
solutions  for  T  >  0. 

The  problem  is  formulated  in  Section  2.  The  numerical  procedure  is 
presented  in  Section  3  and  the  results  are  discussed  in  Section  4  and  5. 

2.  Pormulation 

Let  us  consider  the  steady  two-dimensional  potential  flow  of  an  inviscid 
incompressible  fluid  past  a  bubble  (see  Fig.  1).  The  constant  pressure  in  the 
bubble  is  denoted  by  P^.  We  Introduce  Cartesian  coordinates  with  the  origin  at 
the  top  of  the  bubble  and  we  assume  that  the  bubble  is  symnetric  about  the  x- 
axis.  Gravity  acts  in  the  positive  x-direction. 

We  approximate  the  wake  behind  the  bubble  by  a  free  streamline  model. 
Therefore  the  velocity  is  equal  to  U  on  the  surfaces  SJ  and  S'J*  of  the 
wake.  Inside  the  wake  the  fluid  is  at  rest  and  the  pressure  is  hydrostatic. 

We  define  dimensionless  variables  by  taking  U  as  the  unit  velocity  and  D 
as  the  unit  length.  We  introduce  the  potential  function  fb  and  the  stream 
function  ipb.  The  constant  b  is  chosen  such  that  «  1  at  the  separation 
points  S  and  S'.  Without  loss  of  generality  we  choose  ♦-O  at  x*y“0 
and  ij)  >  0  on  the  surface  of  the  bubble  and  of  the  wake. 

We  denote  the  complex  velocity  by  u  -  iv  and  we  define  the  function 
T  -  i0  by  the  relation 

u  -  iv  *  e^  .  (3) 

We  shall  seek  t  -  10  as  an  analytic  function  of  f  *  ij>  +  il^  in  the  half  plane 
<  0.  The  complex  potential  plane  is  sketched  in  Figure  2.  At  infinity  we 
require  the  velocity  to  be  unity  in  the  x-direction  so  that  t  -  i0  vanishes  at 
infinity  in  view  of  (3). 

On  the  surface  of  the  bubble  the  Bernoulli  equation  and  the  pressure  jump 
due  to  surface  tension  yield 


12  T  b 

—  q  -gx+—  K*B  +  —  on  SAS'  , 
2  p  p 


P  P 
w  b  T 

—  -  —  «  —  K  on  SES'  . 
P  P  P 


Here  q  Is  the  flow  speed,  K  the  curvature  of  the  bubble  surface  counted 
positive  when  the  bubble  is  on  the  concave  side  of  the  surface,  B  the 
Bernoulli  constant  and  p^  the  hydrostatic  pressure  in  the  wake.  In 
dimensionless  variables  (4)  becomes 

e^  -  — ^  X  -  —  e^  0<4  <1  . 


Here  F  and  a  are  the  Froude  and  Weber  number  defined  by  (1)  and  (2) 
respectively . 

On  the  surface  SJ  and  S'J'  of  the  wake,  the  velocity  is  equal  to  U. 

In  dimensionless  variables  this  yields 

T  «  0,  4(  ■  0,  ♦  >  1  .  (7) 

It  is  convenient  to  eliminate  y,  P^  and  B  from  (6)  by  differentiating 


(6)  with  respect  to  4  •  Using  the  relation 

l2L+  ilLm  _J - - 

34,  34,  u  -  iv 


-T-H0 


we  obtain 


e^  e*^ cos  0  -  {e  -  0,  \|(  -  0,  0  <  ♦  <  1  .  (9) 

34  p2  o  3^  3^ 


At  the  top  of  the  bubble,  the  velocity  is  equal  to  zero  and  the  profile  of 
the  free  surface  is  tangent  to  the  x-axis.  This  yields  the  conditions 

T*-«at4»'|>*0  (10) 

0*0  at  4  *4  »  0  .  (11) 

Finally  the  symmetry  of  the  problem  implies 

0-0,  iJ/^O,  4<0.  (12) 

This  completes  the  formulation  of  the  problem  of  determining  the  function 


T  -  i0  and  the  constant  b.  For  ech  value  of  o  and  F,  t  -  10  must  be 
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V  S  ^  S 


analytic  in  the  half-plane  i)>  <  0  and  satisfy  the  boundary  conditions  (7), 
(9)-(12).  This  system  of  equations  is  solved  in  the  next  section.  Once  the 
function  x  -  10  is  known,  the  shape  of  the  upper  suface  SAS'  of  the  bubble 
is  obtained  by  integrating  numerically  (8)  on  tl»  *  0  from  (fr  *  0  to  ^  =  1. 

The  Bernoulli  equation  and  the  continuity  of  the  pressure  across  the 
surface  SJ  and  S'J*  of  the  wake  yield  an  expression  for  inside  the  wake: 

'’w 

—  =B+gx-—  .  (13) 

Upon  substituting  (13)  into  (5)  we  obtain  a  simple  second  order  differential 

equation  for  the  shape  of  the  lower  surface  SES*  of  the  bubble.  Two  boundary 

conditions  for  this  equation  are  obtained  by  imposing  the  continuity  of  the 
tiubbie  profile  at  S  and  S'.  In  the  particular  case  o  =  »  (i.e.  T  =  0) 

lower  surface  of  the  bubble  is  simply  a  straight  line  from  S  to  S’. 


3. 


Numerical  Procedure 
We  define  the  new  variable 

f 


t  by  the  transformation 

1/2  _  _  I)  1. 

■  t'  2i  • 


(14) 


This  c.i ansformation  maps  the  flow  domain  onto  the  inside  of  the  unit  circle. 

The  problem  in  the  complex  t-plane  is  illustrated  in  Fig.  3. 

We  use  the  procedure  derived  by  Vanden-Broeck^  to  investigate  the  effect  of 
surface  tension  on  the  shape  of  fingers  in  a  Hele-Shaw  cell  (see  also  Vanden- 
Broeck^).  Thus  we  define  a  modified  problem  by  replacing  (11)  by 


9=Y,  4,  =  il;  =  0  .  (15) 

Here  y  is  to  be  found  as  part  of  the  solution. 

We  will  solve  the  modified  problem  defined  by  (7),  (9),  (10),  (12)  and  (15) 
for  all  values  of  F  and  a.  Then  we  will  obtain  the  solutions  of  the  original 


problem  by  selecting  among  the  solutions  of  the  modified  problem  those  for  which 


(16) 


Following  Brodetsky'  we  Introduce  the  function  f)(t)  by  the  relation 

T  -  10  *  -  ^  log  -  n(t)  . 

Ti  I  •  t 

The  conditions  (7)  and  (12)  show  that  n<t)  can  be  expressed  as  a  Taylor 
expansion  In  odd  powers  of  t«  Hence 


T  -  10 


I,  V 

n*1 


2n-1 


(17) 


The  function  (17)  satisfies  the  conditions  (7),  (10),  (12)  and  (15).  The 
coefficients  and  b  have  to  be  determined  to  satisfy  the  condition  (9)  on 

the  surface  SAS'  of  the  bubble.  We  use  the  notation  t  *  re^°  so  that  points 


on  SAS*  are  given  by  r“1,  -y<o<^.  Using  (3)  and  (14)  we  have 


dx(q) 

do 


b  sin2oe 


(o) 


cos  6(0) 


(18) 


*  b  sin2oe  ^^*^^8in  0(o)  (19) 

do 

^  M  AV 

where  x(o),  y(o),  T (o )  and  6(0)  denote  the  values  of  x,  y,  t  and  0  on 
the  free  surface  SA. 

We  solve  the  problem  approximately  by  truncating  the  infinite  series  in 
(17)  after  H  -  1  terms.  We  find  H  -  1  coefficients  A^^  and  the  constant 
b  and  y  by  collocation.  Thus  we  introduce  the  N  mesh  points 

Using  (17)-(19)  and  (20)  we  evaluate  t(o^),  0(o^)  and  io  terms  of 

0=0^ 

the  coefficients  A^^  and  the  constants  b  and  y.  Substituting  these 
expressions  into  (9)  we  obtain  N  algebraic  equations  for  the  N  -t-  1  unknowns 
A^,  b  and  y. 

The  last  equation  is  obtained  by  specifying 


(21 ) 


We  solve  this  system  by  Newton  method  for  given  values  of  F  and  a.  Once  this 


system  is  solved,  we  obtain  the  profile  of  the  bubble  by  following  the  procedure 
described  at  the  end  of  Section  2. 

4.  Solutions  Without  Surface  Tension 

We  used  the  scheme  presented  in  Section  3  to  compute  solutions  with  a  >  <» 
(r.e.  T  =  0).  The  numerical  results  were  found  to  be  similar  to  those  obtained 
by  Vanden-Broeck^'®  for  the  flow  past  a  two-dimensional  bubble  in  a  tube.  In 
particular  we  found 

Y  =  j  F  <  ~  0.9  ,  (22) 

Y  =  j  F  =  F^  ~  0.9  ,  (23) 

Y  »  0  F  >  Fj,  ~  0.9  .  (24) 

Ag  F  +  *  the  free  surfaces  collapse  to  the  vertical  line  y  *  0  and  the 

flow  reduces  to  a  uniform  stream.  As  F  0,  the  free  surfaces  approach  the 
horizontal  line  x  =  0  and  the  problem  reduces  to  the  classical  Kirchhoff  model 
for  cavitating  flow  past  a  flat  plate. 

Relations  (22)-(24)  show  that  all  the  solutions  corresponding  to  F  <  F^ 
are  solutions  of  the  original  problem.  The  solutions  for  F  >  F^  are  only 
solutions  of  the  modified  problem.  They  are  characterized  by  a  discontinuity  in 
slope  at  the  apex  of  the  bubble. 

In  Fig.  4  we  present  values  of  the  dimensionless  velocity  at  the 

apex  A  of  the  bubble  vs  F.  The  velocity  is  equal  to  zero  for  F  <  F^ 

and  is  different  from  zero  for  F  >  F^. 

Collins^  performed  experiments  for  two  dimensional  bubbles  rising  in  a 
fluid.  He  observed  that  the  profile  of  the  bubbles  were  close  to  circular 

Q 

arcs.  Using  a  method  due  to  Davies  and  Taylor  ,  he  derived  analytically  an 
approximate  relation  between  the  velocity  0  at  which  the  bubble  rises  and  the 


radius  R  of  the  circular  arc,  namely 


In  dlinenslonless  variables  (25)  becomes 

R  -  4F^  .  (26) 

Collins  found  that  (25)  was  In  good  agreement  with  his  experimental  data.  In 
addition  he  measvired  the  angle  subtended  at  the  center  of  the  arc  of  circle  and 
found  the  approximate  value  105*.  Therefore  the  width  O  of  the  bubble  Is 
related  to  the  radius  R  by  the  simple  formula 

O  -  2  sin  52.5*  R  «  1.6R  .  (27) 

In  dimensionless  variables/  relations  (25)  and  (27)  Imply  the  existence  of 
a  unique  flow  characterised  by 

F  -  ~  0.4  .  (28) 

Our  numerical  results  without  surface  tension  do  not  agree  with  Collins' 
experimental  data:  we  found  a  solution  for  all  values  of  F  <  F^  whereas 
Collins  found  a  unique  solution  for  F  ~  F^.  in  the  next  section  we  show  that 
this  discrepancy  between  theory  and  experiments  Is  removed  by  taking  Into 
account  the  effect  of  surface  tension. 

5.  Solutions  with  Surface  Tension 

We  used  the  numerical  scheme  of  Section  3  to  confute  solutions  of  the 
modified  problem  for  various  values  of  F  and  a. 

In  Fig.  5  we  present  values  of  y  vs  F  for  a  »  40.  as  f  tends  to 
Infinity/  7  approaches  zero.  As  F  approaches  zero/  y  oscillates  often 
around  Fig.  5  suggests  that  there  exists  a  countably  Infinite  number  of 

values  of  F  for  which  ^  solutions  corresponding  to  these  values 

of  F  are  solutions  of  the  original  problem.  As  a  Increases/  the  amplitudes 
and  the  wavelengths  of  the  oscillations  In  Fig.  5  decrease.  These  results  are 
similar  to  those  obtained  in  Ref.  3  for  the  flow  in  a  Hele  Shaw  cell  and  in  Ref. 


4  for  the  flow  past  a  bubble  In  a  tube.  By  analogy  with  those  previous  works  we 


expect  that  the  discrete  set  of  solutions  for  which  Y  ■  will  reduce  to  a 
u:iiq\ie  solution  as  a"  ^  0.  We  shall  denote  by  F*  the  corresponding  value  of 

the  Froude  number* 

In  Ref.  4,  the  value  of  F*  was  obtained  by  following  solutions  with 
7  -  ^  as  a”  ^  I.  (see  Fig.  5  in  Ref.  4).  A  similar  method  was  used  in 

Vief .  3. 

In  the  present  paper,  we  shall  use  a  different  approach  to  determine  the 
value  of  F*.  We  first  recall  that  there  exists  a  branch  of  solutions 
characterized  by  0  <  F  <  F^,  Y  =  ir /2  and  a“^  *  0  (see  Section  4).  On  the 
other  hand  we  expect  to  obtain  a  unique  solution  characterized  by  F  *  F*  and 
Y  =  n /2  as  a  ■►0.  Keller  pointed  out  that  these  two  facts  imply  that 
F  =  F*  is  a  bifurcation  point  on  the  branch  of  0  <  F  <  F^,y  =  if/2,a”^  *  0. 

We  can  therefore  find  the  value  of  F*  by  seeking  numerically  bifurcation 
points  on  the  branch  of  solutions  0<F<F^,Y  *  11/2,  a”'  *  0.  For  this 
purpose  we  use  a  variant  of  the  scheme  of  Section  3  in  which  the  problem  is 
reduced  to  a  system  of  W  +  1  algebraic  equations  for  the  N  +  1  unknowns  b, 
n  =  1,  ...,N  -  1  to  be  solved  for  Y  "  ^  and  a  given  value  of  F.  We 
used  this  scheme  to  compute  the  branch  of  solutions  0  <  F  <  “1^/2, 

=  0.  The  numerical  results  were  found  to  agree  with  those  of  Section  4.  In 
particular  a”  ^  ~  0  for  N  large  (for  example  |®”^|  <  10"’^  for  N  *  40 


•*] 


ind  |a~^|  <  10""^  for  N=60). 


For  each  value  of  F,  we  conputed  the  value  of  the  determinant  of  the 
Jacobian  matrix.  We  found  that  it  vanishes  for  F  ~  0.51.  Therefore  F~  0.51 
is  a  bifurcation  point  on  the  branch  of  solutions  0  <  F  <  F  (Keller^®)  and 


F  ~  0.51. 


The  profile  of  the  bubble  corresponding  to  F  *  F  and  o  =  •  is  shown  in 

*  O 

Fig.  6.  The  broken  line  corresponds  to  a  circular  arc  of  radius  R  =  4(F 
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Pig.  6  shows  that  our  theoretical  profile  is  in  very  good  agreement  with  the 
experimental  result  (26).  However,  the  theoretical  value  F*  ~  0.51  is  larger 
than  the  experimental  value  P^  ~  0.4.  This  is  presumably  due  to  our  simplified 
model  of  the  wake  and  to  the  difficulty  of  realizing  truly  two  dimensional  flows 


experimentally. 
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The  motion  of  a  two-dimensional  bubble  rising  at  a  constant  vel  '.i  v 
an  unbounded  fluid  is  solved  by  series  truncation.  It  is  assumed  thei.  ^ 
wake  of  stagnant  liquid  extending  to  infinity  below  the  bubble.  Both  tl.c 
effects  of  gravity  g  and  surface  tension  T  are  taken  into  account 
shown  that  the  problem  is  characterized  by  a  continuum  of  solutions  i 
and  by  a  discrete  set  of  solutions  when  T  >  0.  In  addition,  a  uniqj 
is  obtained  in  the  limit  as  the  surface  tension  approaches  zero.  Thi 
sponding  profile  of  the  bubble  is  found  to  be  in  good  agreement  with  ivvi. 
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